. V." 





^ « 

r 



MTl-KMSNATION OF THE ELEVATOR POWER 
AND THE DAMPiNG IN PITCH OF THE 
CESSNA I4C AIRPLANE FROM 
fUGHT TESTS 

WILLIAM H. LIVINGSTON 
1950 



rhesis 

.73 




I 



( 



j 





<x» ^ 

a 




€ 



REPORT 



AERONAUTICAL ENGINEERING LABORATORY 



DETERMIIIATIOH OF THE EIEVATOR Pa’.ER AlH) 

TJE DAMPD3G PITCH OF TIE CESSHA UO AIEPLAIE FROM 

FLIGHT TESTS 



by 



PRINCETON UNIVERSITY 
AERONAUTICAL ENGINEERING LABORATORY 



PAGE 2 
REPORT 16 



UUDLEy K 
;maval 

SMOVr* ■' Y 



ACia^^OWLEDGSLENT 

The author wishes to express his appreciation to Professor 
Cotirtland D, Perkins under whose direction this investigation was 
conducted. His ideas, suggestions, and criticisms were most 
helpful throughout the pro j set o 

I(lr. Harry Ashworth, Sr. extended a great deal of assistance 
in preparing the instrumentation, for which the author is indeed 
grateful. 

Special appreciation is due lir. Arthur T. Clemen for his ever 
willing and eager service as observer in conducting the flight tests. 

To all others, whose assistance in so many ways has helped to 
make the project possible, the author's thanks. 



TABLE OF CONTENTS 





Page No 


Title 


1 


Acknowledgement 


2 


Table of Contents 


3 


List of Figures and Illustrations 


4 


Summaary 


5 


Objective 


6 


Introduction 


7 


Description of the Airplane 


9 


Instrumentation 


U 


Weight and Balance 


16 


Resiilts and Discussion 


17 


I • General 


17 


II. The Elevator Power, Cj^^ 


18 


III. The Damping in Pitch, 


26 


Conclusions 


36 


Recommendations 


37 


Preferences 


38 


Curves 


39 


Appendix I; List of Symbols and Abbreviations 


48 


Appendix JI: Airspeed Calibration 


50 



Elevator Position Indicator Calibration 51 
Appendix III; Weight and Balance Calculations 52 
Appendix IV; Photographs 55 



LIST OF FIGURES AMD ILLUSTRATIONS 



Fig. No. Title Page K 

1. Elevator Angle vs Equivalent Airspeed, Forward 39 

Center of Gravity 

2. Elevator Angle vs Equivalent Airspeed, Mid Center 40 

of Gravity 

3. Elevator Angle vs Equivalent Airspeed, Aft Center 41 

of Gravity 

4. Elevator Angle vs Airplane Lift Coefficient, 1-g 42 

curves 

5. Elevator Power vs Airplane Lift Coefficient 43 

6. Elevator Angle vs Airplane Lift Coefficient for 44 

all n's. Forward Center of Gravity 

7. Elevator An^e vs Airplane Lift Coefficient for 45 

all n’s.Mid Center of Gravity 

8. Elevator Angle vs. Airplane Lift Coefficient for 46 

all n's. Aft Center of Gravity 

9. Damping in Pitch vs Airplane Lift Coefficient 47 

II-l. Airspeed Indicator Calibration 50 

II-2. Elevator Position Indicator Calibration 51 

IV-1. The Cessna 140 Airplane 55 

IV-2. Sensitive Accelerometer as Mounted in the Airplane 56 



\ 



SUrJvIARY 



Presented herein are the results of a flight test program 
conducted to detennino two of the longitudinal stability deriva- 
tives of the Cessna 140 airplane. These derivatives are the 
elevator power, , and the damping in pitch, The flight 

test techniques employed to obtain these derivatives are outlined. 
The values of Cjjj^ and thus obtained are compared with the 

values of these parameters computed from the theoretical formulae 
as a means of evaltiating the steady state flight test procedures 
employed. Soixrces of error which arise in the experimental de- 
termination of the derivatives are analyzed and techniques of 
minimizing these errors are indicated. 
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0BJ2CTIVE 



The objective of this investigation has been the determination 
of two longitudinal stability derivatives — the elevator power (Cjq^ ) 
end the damping in pitch of a typical light airplane, the 

Cessna 14C, By comparing the experimental values of these deriva- 
tives derived from flight test data with the values calcu?,ated by 
means of the theoroticcl approach, it in hoped that some evaluation 
ox the flight test teclmlques employed to obtain the derivatives can 
bo ma.de. The resvilts may be of some value to the light aii'plane 
designer in determining the required elevator power for airplanes 
with high damping in pitch. 



DATE im PLACE OF TIE INVSSTIGATIOH 

The investigation was conducted under the direction of the 
Princeton University Aeronautical Engineering Department at Princeton 
Airport, Princeton, New Jersey, from January to May, 1950. 



INTRODTJCTIOts^ 



In recent years the increasing speeds and porfonnance of aircraft 
have necessarily focused lauch attention to the stability and control 
problems Involved in sectoring desirable handling qualities of the airpleine 
in high speed flight. Jn the past, much of the necessary stability and 
control data could be successfully predicted with good accuracy from wind 
tunnel tests. Today, however, the high subsonic, transonic, and super- 
sonic speeds under consideration definitely limit the usefulness of the 
wind tunnel in this field. The difficulties of design and construction 
and the prohibitive costs of such tunnels, coupled with the inaccuracies 
of small scale models, force the aerodynamic ist to turn to other means 
to obtain stability and control data. The other possibility is the air- 
plane itself and this has, in fact, been the modem approach to the 
problem. 

In general, it is necessary to deduce the airplane’s aerodynamic 
characteristics from the response to its controls. I.iuch study is now being 
given to the so-called frequency response of the airplane as a means of 
determining the stability parameters. In this approach the airplane is 
oscillated by an oscillatoiy input to a particular control and its steady 
state frequency response is recorded. In some cases the frequency response 
may be deduced from the transient response to a pulsed input to the control. 
The aerodynamic parameters are then calculated from the frequency response 
data. While these methods will undoubtedly bo of increasing usefulness in 
the future, they are at present somewhat limited by the need for precise 
instrumentation and accurate data reduction techniques. Exact knowledge of 
the airplane's mass and inertia characteristics is required and this in 
itself is no small problem. 
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Y:horo tho flight duration of tho airplane is not lirr.ltod aevoi'oly, 
ac in the csr.o of practically all present day subsonic aii’croft, certain 
fli3iit testing tocliiiquos for obtaining soae of the stability paraioeters 
have been developed tjhich are rtuch cisiplor in theory and. executiori than 
the frequency response methods. In ns.ny coses these techniques are not 
only more accurate but tiieir value is enhanced by loTJer cost, simpler 
instnunentation, and more precise data rodnetion. It is T/ith come of 
these procedvires that this paper will deal. 

The techniques under considerat3.cn involve the deduction of the air- 
plane's elevator power and its damping in pitch from tlie co-cellcd trie 
curves taken in unaccelerated end accelerated flight. As e means of 
evaluating the methods, flight tests of this type were conducted on a 
Cessna 140 airplane, and these two stability derivatives were computed 
from the results. Comparison with results obtained irx computing the de- 
rivatives from the theoretical equations should yield soma insight into 
the accuracy that may be expected in deterrining these der3.vatives from 
flight .tests in this manner. 
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DESCRIPTIOi'] OF THE AIRPL/iNE 



The Cessna 140 is a single engine, high wing, two place, 
personal type monoplane with external bracing and fixed, conventional 
landing gear. The wing is rectangular with rounded tips. It has 
the normal configuration, single vertical tail, Frieze tjrpe ailerons, 
and is equipped with trailing edge plain flaps. 

The airplane is of seni-monocoque all metal construction except 
the wings which are fabric covered. All control surfaces are of 
metal. Figure 17-1 of the appendix illustrates the appearance of 
the airplane as instrumented for these flight tests. 

The following general specifications and dimensions are taken 
from the drawings and reports of the manufacturer except in the case 
whore they have been modified by the installation of special equip- 
ment, In these cases both the old and the new dimension or specifi- 
cation is given. 
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Airplane . general 


# - 


Iilantifacturer 


Cessna Aircraft, Co 


Type 


140 


Recocmended gross weight 
Center of gravity range* 


1450 lbs. 


forwaird limit 


22.8/^ mac. 


aft limit 


30.0^ mac. 


Overall length 


256.50 in. 


Height 

Maximum allowable maneuvering load factor 


74.25 in. 


gross weight 1460 lbs. 


/4.57 to -2.26 


flaps down 40° 

Wing 


/I. 97 to -2.26 


Airfoil section 


lAiCA 24-e2 


Span 


394 iii. 


>? 

Area (total) 


159.29 cq. ft. 


Area (less ailerons) 


145.21 sq. ft. ' 


Aspect ratio 


6.76 


Taper ratio 


1.0 


Chord 

Mean aerodynamic chord 


60.5 in. 


Length 


59.02 in. 


Distance of leading edge back of nose 
reference datxr.: line 


56.58 In. 


Incidence 


/1° 


Dihedral 


/1° 
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Aileron 


1 


Type 


Modified Frieze 


Area 


14. OS sq. ft. 


Span 


74 iJa. 


Chord 
. Travel 


14 in. 


up (fx'Oia uoutral) 


22° 


down (from uouti’al) 


1 >0 


Type 


Plain, Trailing edge 


Ai-ea 


F.736 sq. ft. 


Spaii 


7F.625 in. 


C Lord 


^ 8,0 in. 


Travel (down) 

Horizontal tail trurface 


40° 


Airfoil section 


HACA CCC9 


Ai-ea (including elevators) 


24.35 sq. ft. 


Span 


1C6 in. 


Llaximuin chord 


41*4 in. 


Incidence 


-2.5° 


Dihedral 


0 


Elevator area (total, including tab) 


9.66 sq, ft. 


Elevator span 
Elevator travel 


106 in. 



up (from streauilino with ctabilizer) 20° 

down (from streamlino with stabilizer) 20° 
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Elevator trim tab area 


0 • 5 • ft 0 


Elevator trim tab 'span 


36 in. 


Elevator trim tab mean chord 
Elovator trim tab travel 


5.20 in. 


up (from elevator trailing edge) 


6° • 


down (from elevator trailing edge) 
Vertical tail surface 


330 


Area 


12.420 sq. ft. 


Fin area 


cq. fto 


a pan (to fuselage center line) 


52.2 in. 


Kudder iirea 


5.752 sq. ft. 


Rudder span (maociiiruri) 


49.5 in. 


Rudder travel 

right (fi'OD. streamline with fin) 


16® 


left (from etroaad.ine ’.Tith fin) 
Fuselage 


16® 


Maximum width 


40.0 in. 


Miaximum height 


51.0 in. 


Length (tip of nose to tip of tail) 
Engine (Continental) 


256.50 in. 


Type 


C 85 


Number of Cylinders 


4 


Power Settings Condition 


RPM 


Take off 

liRtP 

Cruise 


2575 

2575 

2400 
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Propeller 

Kanxifactvtrer 

Type 

Diameter 



Sensenich Bros. 
V.'ood, fixed pitch 
79 in. 



AhKONAUIICAL tN(jlNttKIN(j LAbUKAIUKy 



Kth'VJKI 



IlNlSTmiME>:TATION 



The instrumentation required to obtain the necessary data for this 
investigation consisted of an airspeed indicator^ an altimeter, an accelero- 
meter, and an elevator position indicator. The individual instniments 
may be described as follows: 

1. Airspeed Indicator : 

The airspeed was measured with a standard sensitive type airspeed 
indicator connected to a full swiveling pitot-static head. This instru- 
ment was calibrated by means of the speed course method and a calibration 
curve is contained in the appendix. 

2. Altimeter : 

The altitude was measured with a standard sensitive type altimeter. 

3. Accelerometer : 

Since the range of accelerations encoiontered in the tests was small 
and since it was desired to measure increments of acceleration as small 
as 0,25 g, the scale of a standard accelerometer was considered to be 
inadequate. A special accelerometer was constructed and calibrated for 
this investigation. 

This accelerometer consisted of a glass tube approximately 20 inches - 
long, a small colled spring, a steel weight, and an attachment for con- 
necting the spring inside one end of the tube. Two additional weights 
were constructed for use in calibration. These two weights were in the 
ratios of two and three to the weight to be used in the Instrument, 

With the spring motmted in one end of the glass tube, the instrument 
was calibrated by securing the tube in a vertical position, hanging each 
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of the three weights on the spring in tiim and carefully marking 
on the tube the equilibrium position of the weight. In this manner 
the tube was calibrated for accelerations of one, two, and three 
"g's", three *'g's" requiring an extension of the spring of about 
ten inches. Assuming the spring constant to be linear, the fractional 
"g” positions were easily located and inscribed on the tube. 

A photograph of this accelerometer as mounted in the airplane 
may be seen in the Appendix. 

4. Elevator Position Indicator ; 

The elevator position was measured by a 26 volt, 400 cycle a-c 
autosyn type transmitter. This transmitter was linked to the elevator 
control cable as close to the elevator horn as practicable. A cali- 
bration of the autosyn in terms of elevator angle is contained in 
the Appendix. 

All data from the flight tests was visually read and manually recorded 



by the pilot and co-pilot 
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WEIGHT AND BAIjmCE 



All flight tests were made for three center of gravity locations, 
namely, 24.6^, 27.6%, and 30.6% mac. The normal center of gravity 
location for the airplane less pilots is 27.4% for a gross weight of 
1100 pounds. The weight of the two pilots and parachutes was 335 
pounds located at the center of gravity in all three cases. The normal 
operating gross weight for this airplane is 1450 pounds. The airplane 
was ballasted to the desired centers of gravity by means of lead 
weights attached either to the tail wheel arm or to the engine mounts 
as required. As all tests were flown with the fuel tanks approximately 
full, and since all tanks are located near the nonnal center of gravity, 
no corrections to center of gravity position were made to account for 
fuel consumption. 
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RESULTS AND DISCUSSION 



I, GENERAL 

The flight test data required for computing the tT?o desired stability 
derivatives, the elevator power and the damping in pitch, is essentially 
the same in both cases. Determination of the elevator power, requires 

curves of 5^ versus at 1 g for several center of gravity locations. De- 
termination of the damping in pitch, C , requires these same curves for 
at least one center of gravity location (preferably several) and, in addition, 
the curves must be obtained not only at 1 g, but for several values of 
acceleration, i.e., several values of n. 

The data was taken for three center of gravity locations as indicated 
under "freight and Balance”. The speed range over which the tests were run 
extended from 60 mph, to 120 mph, For the 1 g crirves, data was taken 

in steady unaccelerated flight at increments in airspeed of 5 mph. For 
values of n greater than 1.0, the data was taken in accelerated flight in 
steady banked turns at increments of 10 mph in airspeed. 

The airspeed instrument error was negligible j therefore Vj[^ was eq\ial to 
Vq. All flight tests were made at a pressure altitude of 5000 feet and 
since the speed range was low, the effects of compressibility were considered 
to be negligible and V^, was taken equal to V^. The flight test data is 
plotted in Figures 1, 2, and 3 as curves of versus V^. 

In order to eliminate the effects of the elevator tab on required 
elevator deflection, and also to eliminate tab corrections in the calculations, 
the elevator tab was trimmed to zero deflection throughout the tests. 
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The greatest source of error (as will be shown later) in determining 
the desired derivatives, came in tLirough the variation in the power effects. 
Since the Cessna I 40 is equipped with a fixed pitch propeller it was im- 
possible to maintain a constant power setting w hil e varying the spee<|. At 
the higher end of the speed range (100 - 120 mph) it was necessary to 
reduce power rather drastically to avoid over-speeding the engine. The 
RPM was held constant throughout the tests. The cruising RPH of 2400 was 
chosen as the setting and was held at this value for all runs. 

The succeeding sections treat the two derivatives under consideration 
in some detail in attempting to show how the value of the derivative is 
determined from reduction of the flight test data end how this value compares 
in magnitude with the theoretically computed derivative. Sources of error 
are analyzed and finally, an evaluation of the technique for obtaining the 
particular derivative is given, 

II. TIE ELEVATOR POWER . 0^^ 

1, Theory ; 

?i/hen an airplane is in equilibrium the sum of the moments acting on 
the airplane about any chosen reference line must be zero. If, for instance, 
the root leading edge of the wing is selected as the datum, or reference 
line, then the moment of the airplane's weight, W, about this dattim is the 
product of the weight and the distance from the reference line to the center 
of gravity of the airplane (assuming the weight of the airplane to be con- 
centrated at the center of gravity). This moment, then, is given by; 

w • ’‘‘•a. 



(1 
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Now, if a new center of gr*avity position is obtained by shifting 
some portion of the airplane's weight (ballast weight - W^) to a position aft 
of the original center of gravity by a distance, d, the moment of the air- 
plane's weight about the reference line is given by: 



Mref • YV • X 



c,q. 






(^) 



Therefore : 



W • = (w -Wb) + Wb (d + 

W*Xc.(|.^= W* + W^' 

or; 

W (Xc.C^.g- X * W^' d (3) 

Obviously, then, the product of the airplane's weight and the change in 

center of gravity position must equal the applied moment, i.e,: 

w(xc.^.^- Xc,c).J = (4) 

In flight, any moment acting on the airplane in the longitudinal plane 
must be balanced by a moment produced by the deflection of the elevator, . 
That is, the sum of the applied moment and the moment due to elevator de- 
flection must be equal to zero. Therefore, the magnitude of the pitching 
moment coefficient produced by the elevator per degree deflection is a 
measure of the elevator power. The moment produced by the elevator is given 
by: 

A Me * ciH» . ASe (5) 

dSe 

Therefore : 
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At equilibrium the lift must be equal to the weight (L • W). tlaking this 
substitution and putting the equation in coefficient form by dividing by 
qSc : 

Cl • ASe (7) 

0 ^ / 

Solving for : 

C , = _ Cl (Xc.<^- Xc.c^iV 

"'S ASe 

where is now given in per cent of the mean aeix)dynainic chord. 

2. Determination of jij 

Examination of the expression for C^j ^ developed in section 1 above, 
indicates that 0^^ can be computed (as a function of airplane lift coefficient) 
from data obtained in simple flight test procedures. This technique may be 
outlined as follows : 

1. Curves of S* versus Ct are obtained for at least two center of 

* ^a 

gravity locations of the aiiTplane. 

2. For a given Cl determine from these curves the change in elevator 
angle required by the given shift in center of gravity. 

3. Compute Cj,j using equation (8). 

It should be noted that the curves of shown in Figure 4 ar© plotted 
versus the airplane's lift coefficient which is, in general, not equal to the 
lift coefficient of the wing alone. The simple relation derived above for 
determining the elevator power, Cj^g, assumes that Cl^^ ■ That this is 

not necessarily true is easily recognized from the fact that the lift of the 
tail is a definite contribution to the total lift of the airplane and varies 
with elevator deflection. For Cl^ to remain a constant as the center of 
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A 



V 



gravity position is shifted requires that for a change in tail lift due 
to elevator deflection a corresponding increase or decrease in must 
occur. The change in wing angle of attack necessary to produce this lift 
change must in tvim be compensated for by a slight change in elevator de- 
flection. It has been shown, however, that the errors introduced by 
assuming ® normal configurations are not large and probably 

are well within the accuracy of the flight tests. (Reference 3). For this 
reason it was considered justifiable to assume s throughout this 
investigation and, therefore, Cj^gwas calculated using the simplified approach 
as outlined above, 

3* Results : 

A curve of Cj_g versus Cl was computed in the above manner as shown in 
Figure 5. Since three equally spaced center of gravity locations were used 
in the flight tests, the values of 5^ used in computing Cj^gWere taken as 
the average change in elevator angle between the curves of versus Cl for 
the three center of gravity locations. Theoretically, for equal incremental 
changes in center of gravity position these curves will be equally spaced. 

That is, the change in elevator angle required to maintain a given Cl will 
be a constant for even increments of center of gravity shift. Also, theo- 
retical curves of 5^ versus Cl for various center of gravity locations will 
intersect at a common elevator deflection at Cl s 0. It will be noted that 
the experimental curves of Figure 4 do not follow the theory exactly, 
particuleirly when Cl is less than 0,4. The curves obtained for the aft and 
mid center of gravity locations tend to follow the theory closely. However, 



m 
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for lift coefficients below 0,4 the curve for the forward center of 
gravity exhibits a tendency to break down. For this reason the elevator 
power, Cjag , was computed only in the range of airplane lift coefficient 
from 0,40 to 1.20, 

The discrepancy found in the forward center of gravity curve was 
also exhibited in previous flight tests of the Cessna I 40 . (Reference 2) . 
The decrease in required elevator angle at the higher speeds may be due 
to one or more of several possibilities. It has been shown (reference 3) 
that as the center of gravity is moved further and further forward, the 
elevator increment per unit center of gravity change Increases, However, 
the deviation of the curve is rather abrupt below r 0,40 and it is 
believed that the discrepancy may be due in whole, or in part, to certain 
aero-elastic effects not presently accounted for, (Stabilizer deflection 
is among the possibilities.) 



4* Theoretical C alculation of 

For an airplane flying in equilibrium (propeller off) the balance of 
moments (in coefficient form) is given by: 



Cm_ * C 






•f 






(Reference 1,) Differentiating thife expression with respect to 



C9) 



dC„ ^ _ [i^ V d«xt_ 

dSe Vd"*. l-t dSe 

= 'C, = evaVor ef-f ec4-weness 

dSe 
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This expression contains no propeller or power tenns. This omission 
produces no error for the propeller wind-milling condition but will give 
results somewhat low for high power, low speed flight because of the 
increase in dynamic pressure at the tail due to the slipstream. 

Calculation of the theoretical value of the elevator power according 
to this equation is as follows: 

f 

1. ,a. : A+ r - 3.2 

144 X 24.35 

This value for the tail aspect ratio is used to obtain the 
slope of the lift curve of the tail, a^, from Ref. 1, Fig. 5*5. 

a^ r 0.048 per degree 

2, V: V - ^ ^ : 0.428 

c 159.29 X 59.02 



3, Assume s 0.90 



4 . 



r 0.398 

S^. 24.35 

♦ 

Using this value of the ratio of elevator area to horizontal 
tail area the theoretical value of Xt is obtained from Ref. 1, 
Fig. 5-33. 



= 0.59 

Therefore: s (0.048 x 0.428 x 0.90 x 0.59) 

c jjj g — 0 . 0109 
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5. Comi^arison of Theoretical and Experimental Values for : 

The theoretically determined value for the elevator power (Cjn ^ z 
— 0.0109) Is plotted In Figure 5 for purposes of comparison with the 
experimental cuive. The theory asstunes ^ to be a constant based on the 
elimination of power effects. As pointed out above, this assumption leads 
to obvious Inaccuracies in the low speed, high power condition due to the 
increase in the so-called tail efficiency, The tail efficiency is 

defined as the ratio of q at the tail to the free stream q. C^g varies 
directly as and hence increases as Cj^ increases. 

This variation of C^g with Cj^ is quite evident in the experimental 
curve, Cjjjg varying from —0.0135 at Cj^ s 0,4 to —0,0167 at Cj^ s 1.2. 
Comparing these values with the theoretical value for this derivative, the 
experimental value is found to vary from 71.5$ higher than theoretical for 
the low speed, high power condition (Cl » 1.2) to 23.8$ higher than theo- 
retical at cruising speed (Cl = 0,4), 

As pointed out previously, however, with a fixed pitch propeller it 
was impossible to attain constant power settings over the speed range in- ' 
vestigated. In order to avoid over-speeding the engine, power was rather 
sharply reduced at the higher speeds. Such reduction, obviously Induces 
a greater variation in than would normally be expected if the power could 
have been held constant. If this variation in the tail efficiency were 
accurately Icnown and applied in making the theoretical calculation of 0^ g 
the discrepancy between theoretical and experimental values would tend to 
be reduced and probably would approach a constant as the lift coefficient 
increases. varies somewhat linearly with Cj^, Hence, applying this 



UD^ 
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correction to the tail efficiency would not affect the linearity of the 
theoretical cui^e. 

The experimental determination of as performed under the as- 

sumption that s Cp^ also leads to values for the elevator power that 
are somewhat high. However, the use of this assumption yields an error 
that is usually within the accuracy of the flight tests. 

6. Evaluation of Experimental Technique for Obtaining Cjag : 

The stick-fixed trim curves are among the most accurate data that may 
be obtained from flight testing an airplane. These curves are usually 
reduced to plots of elevator angle versus lift coefficient, , for several 
center of gravity locations. As the pilot technique required in these tests 
is simple, the accuracy of the data depends only on how accurately the 
elevator angle and the airspeed can be measured (assuming, of course, that 
the gross weight and the center of gravity location are accurately known). 
Excellent instrumentation is available for meapuring these parameters. 

The greatest errors involved in obtaining the elevator pxjwer with this 
technique come in through the power effects and the variation in For 

an airplane equipped with a controllable pitch propeller, the power can at 
least be held constant throughout the speed range of the investigation, 
nevertheless, it is extremely difficult to accurately predict the variation 
in and herein lies the greatest soiorce of error. The problem is further 
complicated by the fact that for various airplane configurations the hori- 
zontal tail may be positioned such that it is completely in, partially in, 
or completely out of the slipstream. In the latter case, of course, the 
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problem is simplified since is then a constant. For horizontal tails 
only partially submerged in the slipstream it is again difficult to assess 
just what portion of the tail is actually affected by the slipstream. 

In spite of these errors, however, this flight test technique produces 
experin»ntal values of Cjj, g that are in good agreement with theory and ^ 
should prove to be a valuable tool in determining this parameter. 



III. THE DAKPIMG IN PITCH . 

1. Theory ! 

Consider an airplane flying along a straight path in equilibrium at 
some particular lift coefficient, i.e., unaccelerated flight. The balance 
of pitching moments is then given by: 






c.< 



= c 









Oa) 



At each trim point the lift must be equal to the weight and hence 



is a constant. If the airplane's speed is decreased its lift coeffi- 
cient must Increase and the change in pitching moments due to this increase 
must be balanced by the elevator. If the airplane is placed in a steady 
turn at some bank angle the lift must increase in order to provide a vertical 
component of lift equal to the weight. In a turn, however, the lift coef- 
ficient may be thus increased at no decrease in forward speed and the balance 
of pitching moments will be different from those obtained when the lift 
coefficient is increased in level flight. That is, the change in pitching 
moments, that must be balanced by a change in elevator deflection, will be 
different, and therefore the elevator deflection required to produce a given 
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increment in lift coefficient will not be the same as that deflection 
required to produce the same increment in lift coefficient in level flight. 
Obviously, this , difference in elevator deflection is that increment required 
to overcome the damping in pitch of the airplane. 

If the damping derivative is written in the form Cjj, where (d) is the 



"do 



operator 



iWt ) 



( 'C = airplane's time characteristic = 



m 



f>S V 

balance of pitching moments for accelerated flight is given by; 



) , the 



(\Z) 

As shown in the section on elevator power, the elevator deflection 
must balance out any moments acting on the airplane if equilibrivim is to be 
obtained. Having previously determined the elevator power, Cjjj ^ , from the 
1-g trim cxirves, equation (13) offers a method of determining from 
flight test data. 

The airplane lift coefficient is defined as Ct r , At a 

^ /O v2 

given velocity, V, then, the lift coefficient on a 2-g curve is exactly 

twice that on a 1-g curve. If curves of g versus for several values 

of g are taken in steady turns and re plotted as 5^^ versus Cj;^, the increment 

in S between two of these curves (between 1-g and Z-g for example) will be 
o 

that increment in elevator deflection required to overcome the damping in 

pitch of the airplane. If no damping moments existed the ctirves of versus 

Ct for various g's would coincide since a given elevator deflection would 
^a 

produce a given lift coefficient at any acceleration. This assumes, of course, 
that the airplane's static stability does not vary with the different slip- 
stream characteristics encountered at the various power settings used in ob- 
taining the data. 
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If this assumption Is made, these cui^es of S versus Ct for 

e -^a 

various g’s may be reduced to obtain the damping In pitch derivative, 

In the following manner. The airplane's angular velocity In pitch 
may be shown to be: 






(Reference 1.) 

Or, in terms of the operator (d) : 

de = Ci_a ( 1 -j_'\ 

- 2 TV n*; 

For any value of the airplane lift coefficient, then, d© may be easily 
determined. The equilibrium condition is given by: 



Therefore : 



•m<la 



■"'de 



d © + C 






/iSp = 0 



_ _ Cm/; • k.Se. 

de 



( 14 ) 



(IS) 



(ie) 



(»T) 



2, Determination of 

Equation (17) developed above indicates that can be computed 

from flight test data obtained in accelerated flight. The technique may 

be outlined as follows: / 

/ 

1, Curves of versus for sevejfal values of n, including n s 1,0, 

are obtained for one or more center of gravity locations of the 
airplane, 

2, For a particular Cl and a particular n, compute d©, the airplane's 
rate of pitch in terms of the operator (d), 

, At this particular Clj^ determine the L required to accelerate 
the airplane from n z 1,0 to the particular n under consideration. 



3 
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4» Having previously detennined Cj^ g as outlined in the section on 
elevator power, the damping in pitch, is then easily de- 

termined from equatioft (17). 

5. This process is repeated for a number of Ct 's and n*s. Althoxigh, 

^a 

in general, is not a function of lift coefficient, it is con- 

venient to plot versus in order to establish an average 
value for the derivative. The scatter of the computed points is 
usually considerable. However, if a sufficient number of points 
are taken, the trend of the cuinre will be indicated well enough 
that the curve itself can be determined with fairly good accuracy. 
The final value for then. Is the average value of this curve. 



3. Results ; » 

The curves ot S ^ versus V obtained in accelerated flight were replotted 
o y 

as curves of versus Cr as shown in Figures 6, 7, and 8. From these 
e L.Q 

cujrves f'm^g computed for every possible combination of and n as out- 
lined in the preceding section, A plot was then made of versus as 
shown in Figure 9. As expected, considerable scatter of the individual 
points was evident in this plot. However, the trend of the curve is clearly 
indicated. 7?hile the individual points computed are not shown in Figure 9, 
the range of the scatter is indicated by the cross-hatched band. As Gjj^g 
is, in general, not a function of lift coefficient, the experimental curve 
of Figure 9 was considered to be constant over the range of lift coefficient 
investigated. The average magnitude was determined to be — 0,195. 
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4. Theoretical C alculation of 

The total damping in pitch of an airplane is the sum of the damping 
contributions of the various airplane components — fuselage, wing, hori- 
zontal tail, and the propeller. In the normal configured airplane the 
damping due to the tail is by far the largest factor involved. It is the 
usual practice to evaluate the damping in pitch due to the tail and then 
increase the result by some arbitrary factor to take care of all the other 
contributions to the total damping. This arbitrary multiplying factor is 
normally taken equal to 1,10, 

The damping due to the tail occurs as a direct result of the change 
in the effective angle of attack of the tail produced by the pitching 
velocity of the airplane, i,e,, the angular velocity 0, This change in 
is given by: 

A.o<^ = (ie) 

The pitching moment produced by this angle of attack is ; 



Cm = - Ot It V ^ it 



dt y 



Then; 



dc 



m 



d 



= -Qt a-tV it 



\l) V 

Dividing each side by the time characteristic, V : 




-Qt V 

vt 



(\9) 

(za) 



( 2 .') 



•"'de 



= - V ii.J- 

C -U. 



( 22 ) 



From which; 
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To account for the rest of the airplane's damping in pitch multiply 
by 1,10 giving as the final expression: 



Cmds = -1.10 Qt V -(t ^,J_ (^3) 

Calculation of the theoretical value of the^amping in pitch according 



to this equation is as follovJs: 

1, = 0,048 per degree i 

2. V - 0.428 



3. Assume = 0,90 



2.75 per radian 



4. -±L 

C 



" r 

r- 



13,8 X 12 

59.02 



2.8 



m 






S c 



1450 X 12 



(.00205) X 159.29 X 59.02 X 32.2 



28.0 






®de 



28.0 



0.0357 



- - 1.1 X 2.75 X 0.428 x 0.90 x 2.8 x 0.0357 



m. 



‘d© 



= - 0.1165 



5. Conparison of Theoretical and Experimental Values for 

The theoretically determined value for the damping in pitch (Cjn^g - 
- 0.1165) is plotted in Figure 9 to facilitate comparison with the experi- 
mental curve. The value of this parameter as deduced from the flight test 
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data is - 0.195 as shown in Figure 9* The experimental value is 66.8^ 
higher than theoretical. This large discrepancy between experimented and 
theoretical results may be due to one or nrare of several factors. In the 
following paragraphs these possibilities are discussed in some detail. 

In computing the experimental values of at the various Cjj^’s, 
the magnitude of the elevator power, Cjjjg , used in equation 17 to determine 

at a particular was taken directly from the experimental curve of 
Figure 5* As pointed out in the discussion oh elevator power, the experi- 
mental curve of g versus was also somewhat higher than theoretical, ' 
This discrepancy was attributed to the fact that the variation in •It was 
unknown and to the assumption that been shown that if 

this assvimption is not made, that is, the elevator power is computed on the 
basis of the lift coefficient of the wing alone when Cp, is not equal to 
Cp^, then the elevator power will be reduced to a lower value, (Reference 3 ) 
Since varies directly with Cj^ g such a reduction obviously would also 
decrease the magnitude of 

Another important factor that has been assumed negligible in these 
calculations is the variation of the airplane’s static stick-fixed longi- 
tudinal stability with the variation in power. It was assumed that the 
change in elevator angle required to accelerate the airplane to an n greater 
than 1,0 at a partioular Ct , i,e,, the hS between n curves, was entirely 
due to the damping in pitch of the aii*plane. It can be shown that; 



d ^ dCm/ ciCu 
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That is, the elevator angle required to vary the equilibrium lift 
coefficient varies inversely with the elevator power, Cjjj^ , and directly 
with the static stick-fixed longitudinal stability, dCj/dCj^. The as- 
STOiaption, then, that the total between two different n's at a parti- 

cular Cl includes the assumption that the static stability remains a 
constant throughout the test. This is probably a poor assumption in view 
of the fact that in some airplanes there is a wide variation in dC^/dCL 
with variation in power. (This variation is a function of the relative 
positions of the airplane’s thrust line and its center of gravity location 
in the vertical direction). If this is the case, then, a wide variation 
in power as experienced in these flight tests may be expected to Introduce 
considerable error in the value obtained for If the variation of 

the static stability with power is known, or can be accurately predicted, 
such that the change in 5^ required by the static stability change can be 
deducted from the total required to accelerate at a given Cl^> then 

Cnyg can be determined more accurately and its magnitude will be lower for 
the case where static stability increases as power is increased. 

A third source of error which may account for part of the discrepancy 
between theoretical and experimental results lies in the theoretical cal- 

Cr 

culation itself. It will be renumbered that the damping in pitch of the 
airplane was considered to be due primarily to the horizontal tail damping. 
All other damping contributions due to the fuselage, wing, and propeller 
were Iximped together as a multiplying factor of 1,10 to be applied to the 
horizontal tail damping. There is some evidence to indicate that this 
factor of 1,10 may be too low for airplanes having low relative densities. 
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Zimmerinan developed a theoretical formula for in 1935 in which this 
multiplier was taken to be 1»25« (Reference 4.) Since this time the 
relative density of most airplanes has increased considerably and the factor 
of 1.10 is generally considered adequate. For modem light airplanes with 
low wing loadings, however, the relative density factor, yU. , remains low, 
which in Itself tends to increase the damping parameter, since 

varies inversely with jJL . Nevertheless, it is believed that increasing 
the multiplier from 1.10 to 1.25 or higher when working with light airplanes 
may lead to more accurate results in the theoretical calculation of 
Increasing the multiplying factor to 1.25 yields a for the Cessna 140 
of — 0.133. The experimental value obtained from the flight tests is then 
only 46»6/S higher than theoretical as compared with 66.8^5 when the factor 
is 1.10. This new value, of course, is calculated \mder the same assumptions 
as befoire, i.e., the variation in and the change in static stability 
due to power are neglected. 

6. Evaluation of Experimental Technique for Obtaining 

In the preceding sections a flight test technique for evaluating the 
airplane’s damping in pitch has been outlined and discussed in some detail. 

The usual question now arises as to' whether or not the method produces 
results of sufficient accuracy to be practical. 

The pilot technique required to obtain trim curves in accelerated flight 
is relatively simple. The steady banked turn at constant speed is the 
easiest and most precise method of securing this data. The variables involved 
— the airspeed, acceleration, and the elevator deflection — can be measured 
with considerable accuracy. If the elevator power is accurately known, then. 




M 
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it should be possible to reduce the flight test data to a reasonably 
exact value for 

The major problems to be solved in calculating in this manner are 

the variation of with lift coefficient and the variation of static 
stability with power. The sources of error in determining Cjjj g have already 
been analyzed in previous sections. The static stability variation is parti- 
cularly difficult with propeller driven airplanes. If the aircraft is jet 
or rocket powered there may still be i^'portant variations in the stability, 
but usually they are far less serious than for propeller driven airplanes. 

In high speed aircraft, however, there are additional problems arising from 
inaccuracies in the measurement of the extremely small elevator deflections 
required. 

In spite of the difficulties encountered in separating the variables and 
minimizing the errors encountered, the value of this technique for obtaining 
the damping in pitch derivative should not be overlooked. In the past it has 
been necessary to rely on approximations of this parameter. At best, these 
theoretical approximations were in the nature of an "educated guess". Recently, 

more elaborate techniques for obtaining the longitudinal derivatives from the 

/ 

frequency response of the airplane have been developed. can be determined . 

in this manner but the method requires pi^eclsion Instrumentation and data 
reduction. Errors are cumulative throughout the procedure such that by the 
time the value of the damping derivative is reached the results are likely to 
be much less accurate than those obtained from the steady state flight tests. 

It is believed, therefore, that with good pilot technique in obtaining the re- 
quired data, dependable instrumentation, and the elimination of as many of the 
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variables as possible, the steady state flight tests for determining the 
damping in pitch derivative as outlined will prove to be the most accurate 
method of obtaining this derivative for subsonic aircraft. Some refinements, 
such as a means of handling the static stability changes with power, are 
desirable and should be the subjects of futtnre studies. 

CONCLUSIONS 

Pi*om the results of the flight test program conducted on the Cessna 140 
airplane, it is concluded that: 

!• The elevator power, Cjj^ , can be predicted with good accuracy from the 

\ 

unaccelerated trim curves for the airplane provided that; 

a) • Power effects are carefully analyzed in order that errors due to 

variation in the tail efficiency, can be minimized. 

b) . Consideration is given to the fact that the airplane lift coefficient 

may not be zero when the wing lift coefficient goes to zero. The as- 
sumption that CLq s fairly good in most cases but the accuracy 

of the data reduction will be increased somewhat if the assumption 
is not made and the calculations are corrected accordingly. 

2, The damping in pitch, Cjjyg, can be predicted with good accuracy from the 
accelerated trim curves for the airplane provided that; 

a) . The elevator power, Cja g , can be determined accurately. 

b) . The variation in static stick-fixed longitudinal stability with power 

is either small enough to be neglected or can be accurately predicted, 

3. In making the theoi*etical calcuTation of the multiplying factor of 1,10 
used to correct the horizontal tail damping to include the damping of the 
whole airplane may be too low for airplanes having low relative densities. 



% 
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RBCOfjliENDATiaHS 



On the basis of the above conclnslons the following reconmendations 
are made: 

1» The unaccelerated trim curves of the Cessna 140 should be further investi- 
gated. In particular, the 1-g curve for the forward center of gravity 
location should be w-flotm to check it below - 0,4. The behavior 
of this curve should be carefully studied, 

2, A further study of the elevator power of the Cessna 140 should be made 

in which C^^ is computed on the basis of the wing lift coefficient alone 

1 

as outlined in reference 3, 

3* A flight test program should be carried out on the Cessna 140 to determine 
the variation of static stick-fixed longitudinal stability with power in 
order that the computed value of can be corrected for this effect, 

4* A further study should be made in the theoretical determination of 

to investigate the effect of low relative density of the airplane on the 
result. Damping in pitch due to the propeller, wing, and fuselage shoiild 
be studied to determine more accurately what portion of the total damping 
must be attributed to the components. 



AERONAUTICAL ENGINEERING LABORATORY 



REPORT 16C' 



RE FEKEINCES 



1* Perkins and Hage, Airplane Performance . Stability and Control, 

John T/iley and Sons, Inc., Hew York, 1949. 

2* Graham, Dunstan, Longitudinal Stability and Control Flight Tests 
of the Cessna 140 Airplane . Aeronautical Engineering Laboratory, 
Princeton University, Report No. 111. 

3. Perkins, Courtland D., Methods for Obtaining Aerodynamic Data 
Through Steady State Flight Testing . Fart 1« The Longitudinal 

' Derivatives . Aeronautical Engineering Laboratory, Princeton 
University, 1950. 

4. National Advisory Committee for Aeronautics, Technical Report 521, 

An Analysis of Longitudinal Stability in Power-off Flight With 
* 

Charts for Use in Design , by Charles H. Zimmerman. 




# 



1 







t 



I 







Iv 




i 4 ^ 


1 


i'fl 

■' ."T 




;“ 




‘^'i 


i ’ • 





AERONAUTICAL ENGINEERING LABORATORY 



APPENDIX I 

LIST OF SYI.:B0IS' and ABBPEVIATIOt^S 



REPORT 160 



Air density, slugs/cu. ft. 

True airspeed 

Indicated airspeed 

Calibrated airspeed 

Equivalent airspeed 

Elevator deflection, degrees 

Dynamic pressure, Ibs/sq.ft. = 

Acceleration of gravity, 32,2 ft/sec/sec. 

Normal acceleration in units of g, 

Mean aerodynamic choini 

Area 

Weight 

Mass d W/g 

Moment 

Lift 



Airplane lift coefficient i 

Wing lift coefficient ^ 

Slope of horizontal tail lift curve * 

« 

Slope of wing lift curve 
Angle of attack of wing 
Angle of attack of tail 

^t * ^t 

Tail volume coefficient s ;; 

Tail length (airplane c.g, to aerodynamic center of tail' 
Tail efficiency factor s slipstream q/ free stream q 



' <D 
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e 

e 

d 



m 



Elevator effectiveness d o<|./d 
Airplane's time characteristic m 

/»S V 

Airplane's relative density 
Angle of pitch 
Angular rate of pitch 



Operator s 






f 



6 







w 







i 



t 




mi 




I 



4 




D 



I 





I 
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APPENDIX III 

WEIGHT AND BALANCE COMPUTATION 



Scale Calibration 



Scale 


Load in Pounds 


Serial No. 


50 


100 


150 


200 


250 


300 


350 


400 


450 


500 


550 


600 


F 263218 


151 


201.5 


251.5 


301.5 


351.5 


401.5 




502 


552 


602 


652.5 


702.5 


F 263110 


50.5 


100,5 


150.5 


200.5 


251 


301 


351 


401 




501.5 


551.5 


601.5 


F 263111 


44.5 


100 


150 


200 


250 


300 


350.5 


400.5 




501 


551 


601.5 



Weight* and Balance 

Condition: Oil and Fuel Tanks Poll 

Test Equipment Installed 

1. Weight : 



Scale 


W* 


Tare 


Scale 

corr. 


Weight 


Left Wheel 
(F 263218) 


632# 


2.5# 


-102.5# 


527# 


Right Wheel 
(F 263110 


502# 


1.5# 


-1.5# 


499# 


Tail Wheel 
(F 263111) 


127# 


52.5# 


0 


74# 


Gross weight z 1100# 



* Airplane weighed in level flight attitude 




% 
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2* Balance ; 

Datum line (reference) Root leading edge 

Mean aerodynamic chord . 59 •02" 

Root leading edge to mean aerodynamic chord 

leading edge. . . . 0.47" 



Datum line to main gear contact. 3.31" 

Datum line to tail Tfheel contact 201,10" 



a) Normal center of gravity location ; 

^ (Wp^- V-b) Dn-t>Wt Dt (527 4- 499) (3.31) + (74)(201.1) ' 

Gross weight 1100 

x' r 16.64" 

X * x' — Dist, root L.E. to m.a.c. L.E, s 16,64" - 0.47 : 16.17" 
X s ■ 27.4% m.a.c. 

b) Weight of two pilots and equipment z 335 lbs. 

Located at normal center of gravity location. 

?«eight of airplane and crew s 1100 335 z 1435 lbs. 



®) Forward center of gravity location ; 

Added weight: 16 lead weights s 57.0 lbs. 26,5" forward of datum line, 

- (U35 X 16.64) - (57 X 26.5) 

n = 1492 



X, ; ,i4.f99- r..9A7. z 24.6% m.a.c. 

1 59.02 



U.99" 
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d) Mid center of gravity location ; 

Added weight: 1 lead weight z 1*5 lbs., 191" aft of datum line, 

X - (U35 X 16.64) (1.5 x 191.0) _ 

2 " U36.5 



16.80" - 0.47" 

^2 " 59.02 



27.6^ m.a.c. 



e) Aft center of gravity location ; 

Added weight: 4 lead weights r 15.5 lbs., 191" aft of dattim line. 



3 - 1450 



18.55" 



X - 18.55" - 0.47 

3 “ 59.02" 



30 » 6 % m.a.c. 



s 












f 




I' » 



♦ # 





JUL 2 



Bir40tftYv6/ 



181.93 
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